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ABSTRACT 



►^ ' We study relationships between star-formation rate (SFR) and the accretion luminosity and nuclear obscuration of X-ray selected 

r^ , Active Galactic Nuclei (AGNs) using a combination of deep far-infrared (FIR) and X-ray data in three key extragalactic survey fields 

^ . (GOODS-South, GOODS-North and COSMOS), as part of the PACS Evolutionary Probe (PEP) program. The use of three fields 

^— «, ' with differing areas and depths enables us to explore trends between the global FIR luminosity of the AGN hosts and the luminosity 

\f-\ ' of the active nucleus across 4.5 orders of magnitude in AGN luminosity (L^cn) and spanning redshifts from the Local Universe to 

• , z = 2.5. Using imaging from the Herschel/PACS instrument in 2-3 bands, we combine FIR detections and stacks of undetected objects 

CO . to arrive at mean fluxes for subsamples in bins of redshift and X-ray luminosity. We constrain the importance of AGN-heated dust 

^^ ' emission in the FIR and confirm that the majority of the FIR emission of AGNs is produced by cold dust heated by star-formation in 

Cn ' their host galaxies. 

We uncover characteristic trends between the mean FIR luminosity (Leo) and accretion luminosity of AGNs, which depend both 
on L^GN and redshift. At low AGN luminosities, accretion and SFR are uncorrelated at all redshifts, consistent with a scenario 

. , I where most low-luminosity AGNs are primarily fueled by secular processes in their host galaxies. At high AGN luminosities, a 

^^ ' significant correlation is observed between Lgo and Lagn, but only among AGNs at low and moderate redshifts (z < 1). We interpret 

t I I this observation as a sign of the increasing importance of major-mergers in driving both the growth of super-massive black holes 

C^ . (SMBHs) and global star-formation in their hosts at high AGN luminosities. We also find evidence that the enhancement of SFR in 

luminous AGNs weakens or disappears at high redshifts (z > 1) suggesting that the role of mergers is less important at these epochs. 
At all redshifts, we find essentially no relationship between Lgo and nuclear obscuration across five orders of magnitude in obscuring 
Hydrogen column density (A'^), suggesting that various mechanisms are likely to be responsible for obscuring X-rays in active 
galaxies. 

We discuss a broad scenario which can account for these trends: one in which two different modes of AGN fueling operate in the low- 
and high-luminosity regimes of SMBH accretion. We postulate that the dominant mode of accretion among high-luminosity AGNs 
evolves with redshift. Our study, as well as a body of evidence from the literature and emerging knowledge about the properties of 
high redshift galaxies, supports this scenario. 



1. Introduction 
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sive black holes (SMBHs) and their host galaxies. There 
are now well-established and tight coiTelations between the 
properties of nuclear supermassive black holes (SMBHs) 
and the spheroids of thei r host ga laxies (iMagorriari et al.l 



1 998t iFerrarese & Me rtttl 
Marconi &Huntl 120 03 



Sooo 



"looiT 



Ferrarese et alJ 

Ferrarese & Ford 2005; Giiltek in et al.l 

2009t iGebhardt et al.l l2000t iTremaineet al.ll2002D . The strong 



evolution in the cosmic sta r-formation rate (SFR) density (e.g., 
iHopkins & BeacomI l2006h is also reflected in the evolution 
of the number density of bright Quasi-stellar Objects (QSOs) 
and the cosmic SMBH accretion rate density since z ~ 3 
(Bovle&Terlevich'1998"; 'Ueda et al."2003'; Croom et al. '2004'; 
[Silverman et al. 2008: Aird et a l. 2010; Assefetal. 20 11). In 
addition, AGNs, like star- forming galaxies (Fontan otet alj 
|2009 (). display a form of 'down-sizing', by which the population 
that dominates either nuclear or star-forming activity moves to 
systems at lower luminosities or masse s towards later epochs 



dCowie et all 120031: iFiore et al.1 120031: iHasinger et aTl 120051; 



20031: ^ 
11120071^ 



iBongiorno et al.ll200 

Observational studies of the relationship between host star- 
formation and AGN activity are an essential test of the de- 
gree and importance of galaxy-SMBH co-evolution scenarios. 
Traditionally, such studies of star-formation in AGNs have con- 
centrated on interesting subsamples, usually luminous or rare 
types, such as bright QSOs, radio-loud AGNs or bright local 
Seyfeit galaxies. Many of these types of AG Ns ar e associated 



with o n-going 



1998; Omont et al. 



starb ursts (IRowan-Robinson 
1^031; Prid dev et all 



Haas et al 



-- 1995 

2003^, IJahn keetal. 



2004 IShi et alj|2005l; 1schweitzer et al.ll2006l:lN"etzer et al. ll200 



galax y (ULIRG) regirne (e.g ., |Pi 

120051; [Alexander etai] 120081; ICoppin et al.1 12008). 



Greene et al.l 12009). In general, one finds clear star- formation 
signatures in a sizable fraction of bright QSOs (as much 
as 80%, depending on selection), with some reaching infra- 
red (IR) luminosities placing them in the ultra-lumin ous IR 

Page et al. 2004; Steve ns et al] 
In a sim- 
ilar vein, a significant fraction of ULIRGs/sub-mm galax- 
ies show clear AGN signatures, though not always at the 
level of QSOs ('Sanders et al.'IgS?, 'Sanders &Mirabel|ll996l 
iGenzel e t al. 1998; Canalizo & Stockton 2001). These sfiidies 
support the notion of a connection between galaxy mergers and 
black hole growth, but, by themselves, they do not explore star- 
formation in the bulk of AGNs. 

With the advent of large systematic galaxy surveys, several 
studies have turned to complete, well-defined samples of AGNs 
selected by multiple means and covering a large dynamic range 
in AGN luminosity and obscuration. Star formation in these 
AGN hosts is revealed typically through photo metric or spectro- 
scopic signatures at optical or UV wavelengths. Kauffma nn et alj 
(|2003b calibrated the star-formation history of emission-line se- 
lected AGNs in the local Universe using spectra from Sloan 
Digital Sky Survey (SDSS) and found that the mean SFR and 
typical ages of AGN host galaxies are respectively higher and 
younger than those of inactive galaxies. Post-starburst spectro- 
scopic signatures are also found to be strong in AGN hosts 
dWild et al.ll2007h . while there is evidence in some AGNs that 
nuclear ac tivity follow s a few 100 M yr after a strong starburst 
(iDavies et al. 2007; W ild et alJl20iol) . 

At hi gher redshifts, multi-w avelength sur veys such as 
GOODS dDickinson et alJ l2003l). CO SMOS dScoviUe et alJ 



2007h and A E GIS dDavis et alj feoOTh played a major role 



Silverman et al.l(l2009h . using the SF sensitive [O II] A3727 emis- 
sion line, found wide-spread SF in low and moderate lumi- 
nosity AGN hosts at levels comparable to inactive galaxies of 
similar stellar mass. AGN host colors paint a similar picture 



dXue et al.ll2010l : lRosario et al.l201 ih . Clustering studies suggest 
that AGNs are associated more with red or evolved galaxies, 
rather than the bulk of star-forming galaxies (Coil et al. 200^. 
These studies generally imply that low-luminosity AGNs are 
hosted by a special category of galaxy that is quite evolved 
but still formi ng stars. Since SFR correlates quite strongly with 
steflar mass (Noeske et al.'2007':'Dad di et al.ll2007t ISantini et al.1 
20 09) and AGN hosts are p referentially hosted by massive galax- 
ies dAlonso-Herrero et al.l 12008). stellar mass selection effects 
can bias any infe rences abou t SF in AGN h osts (ISilverman et al.l 
2009; Xu eetal.l l2010: Rosa rio etal.ll2011l) . When taken into ac- 
count, active and inactive galaxies are found to be quite similar 
in their average SF properties. 

A caveat of most studies of star-formation in AGN hosts is 
that they tend to concentrate on low and moderate luminosity 
AGNs. This is because the emission of the active nucleus can 
be a strong contaminant of SF tracers at optical and UV wave- 
lengths. Together with the difficulties stemming from differen- 
tial dust obscuration and reddening between nuclear and galaxy 
emission, this makes it next to impossible to accurately deter- 
mine SF properties of luminous AGNs using these data. Tracers 
in the IR are useful since extinction is rarely a major issue at 
these wavelengths (except around the 10 //m Si absorption fea- 
ture). The IR continuum comes primarily from reprocessed short 
wavelength radiation and the total IR luminosity (generally in- 
tegrated over the range of 8-1000 jum) is considered to be a 
good calorimeter of the total bolometric output of both AGNs 
and SF systems. Since the mid-infrared (MIR) has traditionally 
been more accessible to observation from early generations of 
ground and space-based instruments, several studies of SF in 
AGNs have been based on the measurement of the prominent 
PAH emis sion features that are excited in star-forming envi - 
ronments (iNetzer et al.ll2007t iLutz et al.ll2008t IShi et al.ll2009l) . 
QSOs frequently show stronger PAH features than typical inac- 
tive galaxies, though evolution in the properties of the massive 
galaxy population was generally not taken into account in these 
studies. Another limitation is that MIR spectroscopy, on which 
PAH measurements depend, is only sensitive enough for rather 
IR-bright galaxies, which are not representative of the broad 
population of AGN hosts. 

Wide are a galaxy surveys fr om the new Herschel Space 
Observatory (IPilbratt et al.l 1201 Ol) have opened up the ability 
to carefully measure the rest-frame far infrared (FIR) contin- 
uum luminosity of distant galaxies and AGNs. Several stud- 
ies show that the FIR continuum is the best and least con- 
taminated tracer of SFR in massive galaxies, which are dusty 
enough for most SF luminosity to be reprocessed into the FIR. 
In a first study of X-ray selected AGNs in the GOODS-N field, 
IShao et a l. (2010) found that the mean relationship between the 
Lacn and global host galaxy SF is a function of the accretion 
luminosity, such that luminous AGNs show a SF-Lacn connec- 
tion that is ab sent in lowe r luminosity systems. Using a similar 
approach, Mullanevet al. (2012), using deeper FIR imaging in 
both GOODS fields, find that the mean SFRs of X-ray sources 
is indistinguishable from those of typical inactive star-forming 
galaxies, though the study only concentrated on low to moderate 
luminosity AGNs. Alternatively, Hatziminaoglou et al. (20iy) 
find that FIR luminosity is correlated with AGN luminosity over 
many orders of magnitude in both quantities, at apparent odds 
with aspects of both previous results. 

Many of the differences between these works can be recon- 
ciled by studying the FIR properties of a broad sample of AGNs 
across many epochs. We present the largest study to-date of the 
mean SF properties of AGNs, in terms of the range of accretion 
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Fig. 1. Rest-frame 2-10 keV (hard band) X-ray luminosity Lx vs. hydrogen obscuring column Nh for X-ray sources in this study. 
The upper four panels shows sources in COSMOS, while the lower panels show sources from both GOODS fields combined. Each 
panel is for a different bin in redshift, increasing left to right. In the three higher redshift bins, the location of the X-ray sources from 
the lowest (0.2 < z < 0.5) bin are shown for reference as small red points, to clearly bring out the redshift trends in the sample. 
Sources with Nh consistent with only Galactic obscuration are shown as left-pointing arrows. The fiducial limit of Lx= 10"*^ erg 
s ' used to define a clean AGN sample (employed in §4.3) is shown by a dashed line. The nominal Compton thick regime of X-ray 
obscuration (Nh > 10^"* cm"^) is shown as a shaded region in each panel. 



luminosity (5.5 orders of magnitude), SFR and redshift (Local to 
z - 2.5). This is the second paper of a series that combines deep 
FIR images and the best current X-ray analysis in multiple deep 
fields to determine the SF characteristics of a large sample of 
AGNs spanning the redshifts 0.2 < z < 2.5. Using techniques of 
FIR stacking, we have constrained the mean FIR luminosities for 
active galaxies across redshift, AGN luminosity, obscuration and 
host mass. In the first paper of the series (.Santini et al..2012.) . we 
studied the mean SFR of AGN hosts as function of stellar mass 
and redshift, finding a slight enhancement compared to inactive 
galaxies with the same stellar mass distribution. In th is work, 
we build on early work from our team (IShao etal.ll2010l) and un- 
cover trends with AGN luminosity and obscuration. The paper is 
organised as follows: §2 describes the datasets used and outlines 
our methods (which are detailed in the Appendix), §3 deals with 
AGN emission in the FIR as a source of contamination, §4 out- 
lines our main trends, fits and tests for biases. In §5, we discuss 
our results in the context of galaxy evolution. 

We assume a standard A-CDM Concordance cosmology, 
with Hq = 70 km s"' Mpc"' . 

2. Datasets, Sample Selection and Methods 

This work exploits deep X-ray and FIR imaging, as well as ex- 
tensive ancillary data, in three well-studied extragalactic fields: 
GOODS-South, GOODS-North and COSMOS. The depth of the 
two GOODS fields in the X-ray, FIR and several other bands 
is essential to probe faint and high-redshift AGNs. The shal- 
lower, wide-area COSMOS field provides us with good statistics 
among the more luminous sources, which are rare in the smaller 
GOODS fields. 



2.1. X-ray catalogs 

Detailed X-ray point source catalogs were compiled for all three 
fields. These provide rest-frame X-ray luminosities and obscur- 
ing columns, source classification, optical associations and red- 
shifts. Where possible, spectroscopic redshifts were used. In the 
absence of these, photometric redshifts specifically suited for 
AGN hosts, based on optimised AGN or hybrid templates, were 
compiled. 

For GOODS-S we used the iLuo et al.l (12008 ') catalog from 
the 2 Msec Chandra Deep Field-South imaging program. A 
small set of n ew redshifts from recent spectroscopic efforts 
(ISilverman et al. 20 10) were incl uded, and photom etric redshifts 
were compiled from the work of lLuo et alT(l2010l) . 

For GOODS-N we used the 2 Msec Chandra Deep Field- 
North catalog of Alexander et al. (2003) and an updated version 
of the classifi cation into AGN a nd other t ypes of X-ray so urces 
compiled bv .Bauer et al.l(l2004 (see also .Shaoet al.ll2010l) . 

Finally, for COSMOS, we used the XMM-Newton catalog 
compiled by Cappelluti et alJ (l2009h and the X-ray optical asso- 
ciations and derived properties presented in iBrusa et aP (l2010l) 
in its updated version, which includes new redshifts from on- 
going (mostly DEIMOS/KECK) spectroscopic campaigns and a 
few changes in the published redshifts and /or spectroscopic clas - 
sifications from a re-analysis of the data (Mainieri et al." 201l1) . 
Photometric redshifts are as derived by Salvato et al..(.2009h . 

Rest-frame X-ray luminosities in the hard band (2-10 keV) 
were estimated for all our AGNs, using spectral fits for X- 
ray sources with sufficient counts or using scalings based on 
hardness ratios for faint X-ray sources. Where possible, the 
Hydrogen column densities of obscuring gas towards the X-ray 
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sources (Nh) were also estimated. The X-ray luminosities were 
then corrected to account for the level of obscuration, yielding 
absorption-corrected hard-band luminosities (Lx hereafter). It 
must be kept in mind that Nh estimates, especially for obscured 
AGNs at higher redshifts, can be fairly uncertain. X-ray lumi- 
nosities in the hard band, however, are generally more secure, 
though, of course, for heavily obscured sources, the uncertainty 
in Nh will translate into a greater uncertainty in Lx- 

In this work, we consider the properties of all X-ray sources 
with Lx greater than 10^' erg s"'. At the low-luminosity end 
(Lx< 10^^ erg s~'), there may be some contamination from 
very luminous star-bursts where soft-band emission from star- 
formation related processes (SNe, high-mass X-ray binaries) can 
contaminate the sample. Since we study trends with AGN lu- 
minosity in our main analysis, we can easily test for the im- 
portance of such starburst contamination in our sample. In gen- 
eral, we do not find that low luminosity X-ray sources are more 
strongly star-forming than moderate luminosity sources, except 
in the redshift range 0.2 - 0.5. We account for such effects in our 
analysis. 

A sense of the differing depths of the X-ray data in the three 
fields may be gained from Fig. [T] which plots Lx against the 
obscuring column Nh for the X-ray sources from the three cat- 
alogs used in this analysis. We split the sample into a set of 4 
redshift bins: 0.2 < z < 0.5, 0.5 < z < 0.8,0.8 < z < 1.5 and 

1 .5 < z < 2.5. The fields are complementary in terms of their X- 
ray properties: the deep GOODS fields allow a good sampling of 
faint AGN, while in the large-area COSMOS field we probe the 
rarer and more powerful sources. 

The X-ray sources also display strong trends of Nh with red- 
shift, especially in th e COSMOS field. This is governed mostly 
by selection effects dTozzi et al.l 2006). At all redshifts, a sub- 
stantial fraction of AGNs have absorbing columns consistent 
with no intrinsic absorption (i.e, only absorption by gas in our 
own Galaxy). Towards higher redshifts, the rest-frame soft X-ray 
energies redshift out of the detection bands of both Chandra and 
XMM-Newton, leading to greatly reduced accuracy in the mea- 
surement of Nh in weakly absorbed sources, as well as smaller 
range of estimated Nh values from most spectral fitting meth- 
ods. This explains the large number of sources clustered around 
log Nh = 20.5 cm"^ in the upper panels. In addition, absorbed 
low-luminosity AGNs are badly incomplete in flux-limited X- 
ray surveys, leading to an apparent correlation between Nh and 
Lx, which can be seen most clearly in the high redshift panels 
of Fig. [1] While the span of Lx in any given redshift bin is fairly 
broad, the span of Nh that we can probe in a bin can be quite 
limited. 

2.2. PACS Imaging and catalogs 

Our FIR data come from observations by the PACS instrument 
dPog litsch et al. 2010) on board the Herschel Sp ace Observatory , 
as part of the PACS Evolutionary Probe (PEP Lutz et al."2011) 
survey. Observations of GOODS-S were cariied out in all three 
PACS bands (70, 100 and 160 jum), while the other two fields 
were only observed in the two long wavelength bands. We made 
use of PACS catalogs extracted using prior knowledge of the 
positions and fluxes of sources detected in deep archival IRAC 

3.6 /im and MIPS 24 //m imaging in these fields. This allows 
us to deblend PACS sources in images characterized by a large 
PSF, especially in crowded fields, and improve the complete- 
ness of faint sources at the detection limit. Fluxes in GOODS-S 
reach 1.1, 1.2 and 2.4 mJy at 3cr in the 70, 100 and 160 fim 
bands respectively. In the other two fields the 3cr flux limit at 



100/160 ju m is 3.0/5.7 mJy i n GOODS-N and 5.0/10.2 mJy in 
COSMOS (lBertaetal.ll201lh . Detailed information on the PEP 
survey, observe d fields, da t a proc essing and source extraction 
may be found in lLutz et alJ (1201 ih . 

X-ray sources were matched to PACS sources through their 
optical counterparts. In most cases, a robust crossmatch could 
be made between the optical counterpart of an X-ray source and 
a 24 fim source in the catalogs used to extract PACS fluxes. In 
the few cases where the crossmatches were not robust, either due 
to uncertain associations between optical and 24 fim sources, or 
due to multiplicity in a match, we visually examined the matches 
on the 24 fim images to verify their quality and discard unreliable 
matches. The technique used to associate PACS fluxes with op- 
tical counterparts could in principle introduce some biases in the 
determination of the mean FIR luminosity, since AGNs are typi- 
cally bri ghter than inactive g alaxies at MIR wavelengths (e.g, 24 
fim). In lSantini et al.l (l2012h . we show that this bias is unimpor- 
tant. 



2.3. Swift BAT sample 

For a comparison to local X-ray selected AGNs, we also com- 
piled an unbiased sample of extremely hard X-ray-selected AGN 
(15-1 50 keV band) from the 39 month Palermo Swift-BAT cat- 
alog (Cusumano et al.'2009). This sample has also been used in 
Lutz et al. (2010) and Shao et al. (2010). 

From the Swift-BAT catalog, we selected sources classified 
as Seyferts, LINERs, quasars, and other AGNs. We excluded ob- 
jects with evidence for a strong non-thermal contribution to the 
far-infrared (on the basis of the NED SED) or those classified 
as blazars. In addition, objects at Galactic latitude \b\ < 15 and 
objects at redshift z > 0.3 were also removed. For the remain- 
ing 293 AGNs we used the IRAS Faint Source Catalog 60 fim 
detections where available, otherwise we used SCANPf] to ob- 
tain 60 fim measurements for faint or individually nondetected 
objects. We calculated rest-frame 2-10 keV luminosities extrap- 
olating from the BAT fluxes and the redshift, assuming an AGN 
photon index of 1.8. The sources were averaged in 7 bins of Lx, 
with sufficient statistics in each bin. 



2.4. Methods 

AGNs in all three deep fields were divided into four intervals in 
redshift (0.2-0.5, 0.5-0.8, 0.8-1.5, 1.5-2.5), spanning the epochs 
where most black hole growth in the Universe has occurred. In 
each redshift interval, the AGNs were further divided by Lx 
into five bins (41-42, 42-43, 43-44, 44-45, 45-45.5 in logL^)- 
An estimate of the mean FIR luminosity (vLy at 60 fim) in 
each bin was made by combining PACS photometry for FIR de- 
tected AGNs and PACS stacks for FIR-undetected AGNs in each 
field, and then appropriately combining measurements from all 
fields together for a single measurement per bin (L^o)- Details 
of the PACS photometry, stacking, error estimation and methods 
of combination between and across fields are developed in the 
Appendix. 

3. AGN emission in thie FIR 

The key results in this work are based on the notion that rest- 
frame FIR emission is a proxy of SF activity. A number of pre- 
vious studies suppor t our assumption based on different grounds . 
Many authors (e.g. ISchweitzer et al.1 l2006t iNetzer et all 120071: 



http://scanpiops.ipac.caltech.edu/applications/Scanpi/index.html 
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iLutz et ani2008l) showed a strong correlation between FIR lu- 
minosity and SFR tracers, such as PAH emission features, both 
in local and high redshift bright {Lacn I Lfir up to ~ 10) QSOs. 
In this section, we explore the degree to which AGN emission 
may contribute at FIR wavelengths, in particular at 60 //m. 



3.1. Predictions from Models of AGN-tieated dust 

The nature of emission from AGN-he ated dust has 
the topic of several theoretical studies ("Pier & Krolik| 



been 



Granato & Danese 1994; Efstathiou & Rowan-Robinson 
isfenkova et al.ll2002l: iFritz et al.ll2006l: ISchartmann et all 



1992; 



19951 
200811 



Most models assume that the dust is distributed in a 'torus', i.e, 
a roughly axisymmetric structure with an finite inner radius, de- 
termined by the sublimation temperature of the grains. Smooth 
and clumpy torus geometries models with a range of parameters 
have been studied in the literature. 

A consistent feature of most torus models is a broad peak 
in the continuum emission at MIR wavelengths with a sharp 
dropofFbeyond a certain break wavelength. The MIR SED shape 
shows considerable variation in different models, due to silicate 
emission and absorption features at 10 fim which depend on the 
geometry of the dusty torus and its optical depth in the MIR. The 
typical FIR SED is smoother and more uniform among most 
models. However, a small subset of models, characterised by 
high MIR optical depths (or dust columns) and large torus sizes, 
peak at fairly long wavelengths, approaching 60 /im jFritz et al.l 
12006 !). In such models, the intrinsic IR AGN SED is quite 'cold' 
and a substantial fraction of the reprocessed AGN power can 
come out in the FIR. 

In addition to torus emission, AGN light can also be repro- 
cessed by dust on kpc scales, in the Narrow Line Region (NLR). 
NLR dust is typically cooler than torus dust and has a longer 
tail towards FIR wavelengths (Schweitzer et al. 2008; Mor et al. 
[2009 ). The relative strength of this component with respect to 
the torus emission is a function of the coveri ng factor of NLR 
clouds, which tends to be quite low (~ 0.1; iMor et al.l 120091) . 
indicating that the NLR contribution to the total AGN IR lumi- 
nosity is likely to be small. 

The theoretical approach of modelers is useful to understand 
the possible range of intrinsic AGN SED shapes. However, the 
'torus' and extended circumnuclear environment of the SMBH, 
which governs the intrinsic SED, has only been studied in a few 
nearby AGNs, due to the limited spatial resolution of MIR and 
FIR imagers, and little is known of true variations in their geom- 
etry and phase structure. 

Other works have taken an empirical approach towards cali- 
brating intrinsic AGN SEDs. Such analyses rely on some knowl- 
edge of the range of SEDs of star-forming galaxies, which are 
better constrained from studies of inactive systems. By subtract- 
ing IR templates of star-forming systems from SEDs of active 
galaxies, one may reveal the residual excess whic h must arise 
from A GN-heated dust. Adopting this approach, iNetzer et al.l 
(l2007h . in a study of PG QSOs with bolometric luminosities 
Lagn in the range of lO^'^-lO'**'^ erg s"', derive a mean intrin- 
sic AGN SED that shows strong Si emission features and drops 
rapidly beyond ~ 20 fim. The shape of the mean SED does not 
appear to depend on the level of SF, but does shows considerable 
intrinsic variation from object to object, especially at NIR-MIR 
wavelengths. 

This work was extended recently by 'Mulla nev et al.l (1201 ih . 
which explores intrinsic AGN SEDs in a sample of 1 1 local low- 
to-moderate luminosity hard X-ray selected AGNs with low lev- 
els of SF. Using a set of 5 templates which encompass most of 



the range of star-forming galaxies, they isolate and model the 
AGN emission using a smooth broken power-law model, with 
a sharp long-wavelength dropoff model ed as a modified black- 
body. The empirical AGN SEDs from iMullanev et aH (1201 ih 
show significant intrinsic variation both in MIR and FIR, with 
some peaking as far red as 50 /im. The study also uncovers a de- 
pendence of SED shape on AGN luminosity, such that more lu- 
minous AGNs have, on average, a steeper long-wavelength drop- 
off than lower luminosity A GNs. The long-wavele ngth slope of 
the mean AGN SED from Mullanev et al.l (l20Ilh is shall ower 
than that of the mean QSO SED from lNetzer et al.l (l2007b . but 
this c ould be due to t he low er characteristic AGN luminosity 
of the iMuUanev et alJ (l201ll) sample. Having said this, detailed 
SED modeling of a sample of Type I AGNs, covering a wide 
range in AGN luminosity and including contributions from NLR 
dust emission and hot dust from the outer Broad-Li ne Region, 
seems to suggest that the steep FIR dropoff in the Netz er et al] 
(2007.) SED appli es even to lower lum inosity systems and may 
be more universal (iMor & Netzerl2012l) . Given the recent nature 
of the latter paper, we choose in this paper, to use the SEDs of 
Mullanev et al. (2011) and Netzer et al. (2007), but caution the 
reader that this suite of SEDs may overestimate the contribution 
of t he AGN at FIR wavele ngths . 

iMullanev et al.l (1201 lb also compared their empirical AGN 
SEDs with those predicted from torus models. The empirical 
SEDs shapes show a much smaller sca tter than those of the 
smooth torus models of iFritz et al.l (l2006h . but are more compa- 
ra bly matched to the expect ed range of the clumpy torus models 
of ISchartmann et aP (l2008h . This may arise simply from differ- 
ences in the range of parameter space explored by the modelers - 
the latter clumpy models are dependent on expensive 3-D model- 
ing and e xplore a smaller ra nge of torus sizes and optical depths, 
while the lFritz et al.l (l2006l) models with the largest variation are 
those with large sizes and MIR optical depths. 

3.2. FIR colors of AGNs and inactive galaxies 

A simple approach towards exploring the importance of AGN 
contamination at FIR wavelengths is to compare the FIR colors 
of AGNs against the general population of galaxies. The signif- 
icantly warmer temperatures of AGN-heated dust gives flat or 
falling SEDs at long IR wavelengths, which leads to bluer FIR 
colors among galaxies with a substantial AGN component in the 
FIR. 

In the top panels of Fig. |2l we compare IR colors of PACS- 
detected X-ray AGNs and all PACS-detected galaxies in the 
COSMOS field in four redshift bins. No explicit selection has 
been applied to the PACS-detected galaxies in this plot, except 
that they have a well-defined spectroscopic or photometric red- 
shift. The AGN points are colored by their hard-band X-ray lu- 
minosity (log Lx) using a rainbow color stretch from black (log 
Lx- 42.0 erg s"') to red (log Lx= 45.5 erg s"'), more than 3.5 
orders of magnitude in AGN accretion luminosity. The low end 
of the Lx range is a luminosity threshold which is popularly used 
to identify pure AGN samples, while the upper end overlaps the 
luminosity range of moderately bright QSOs. 

The X-axis of Fig.|2]is the 24 fim to 160 //m color (m24-i6o), 
a MIR to FIR flux ratio, sensitive to warm AGN dust emission. 
At z > 1, the MIPS 24 fim band is also influenced by silicate 
absorption, while, at z > 2, PAH emission contributes as well. 
The FIR 100 /im to 160 fim color on the Y-axis (mioo-ieo) mea- 
sures the FIR continuum slope, down to the rest-frame 30 fim 
at z = 2.5. The flux measurements in these two PACS bands 
are used to derive the rest-frame 60 fim luminosity in all four 
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Fig. 2. Observed IR colors of X-ray AGNs plotted in bins of redshift. A MIR flux ratio (the flux at 24//m to the flux at 160 /im, 
expressed as a magnitude) is plotted against an FIR flux ratio (the flux at 100 //m to the flux at 160 yum, also as a magnitude). The 
shaded grey regions mark out parts of the diagram occupied by star-forming galaxies with an AGN fraction of > 50% at rest-frame 
60 yL(m (§3.2). Top row panels: PACS-detected AGNs from the COSMOS field are shown as star-shaped points, with colors that 
represent their hard-band X-ray luminosity Lx (shown in the color bar at top). The small black dots are all galaxies that are detected 
in the COSMOS PACS maps in the same redshift bins. Objects with colors right and below the dashed line are unlikely to be detected 
because of the photometric limits of the MIR and FIR catalogs. Bottom row panels: Stacked measurements for PACS -undetected 
AGNs from all three fields are shown as star-shaped points, colored by Lx- For reference, the large black points show the locations 
of PACS-detected COSMOS AGNs (the same as the colored points in the top panels). PACS-undetected AGNs have, on average, 
warmer IR colors and may be more AGN dominated in the FIR. 



redshift bins. Therefore, understanding the level of contamina- 
tion from AGN light in these bands is particularly relevant to the 
later discussion on FIR trends with AGN luminosity. Both colors 
are in units of magnitudes; therefore, smaller values (bluer col- 
ors) imply systematically warmer SEDs. We will use the term 
'warmer' to refer to bluer IR colors to avoid any confusion with 
the colors used in the plotting symbols. 

Since the PACS bands probe diff'erent rest-frame wave- 
lengths with redshift, we break the sample into our four standard 
redshift bins to reduce the eff'ect of k-corrections. 

From the Figure, one may notice that the AGNs on aver- 
age show warmer m24-i60 than normal (i.e. X-ray undetected) 
galaxies with FIR detections, though the bulk of AGNs at any 
given redshift have colors that are consistent with the field pop- 
ulation (see a l so the discussion of mid-IR spectra in section 8 of 
iNordon et al.l (I20I2 ')). The most luminous AGNs typically have 
the warmest 24-100 colors. 

The mioo-160 distributions of the AGNs and normal galax- 
ies are m ore similar than their m24-i6o distributions (also see 
iHatzimin aoglou et al.ll2OI0l) . but some small differences are no- 
ticeable. In Fig. [3]we compare the distributions of mioo-ieo be- 
tween AGNs and normal galaxies. At face value, we expect that 
the largest differences between the two sets of histograms should 
be at high redshift. This is because the PACS 100 and 160 ycca. 



bands trace increasingly shorter rest-frame wavelengths, mak- 
ing the AGN emission more prominent with redshift. However, 
Fig. |3] shows that this is not the case. Quite surprisingly, the 
biggest difference between the distributions of AGNs and normal 
galaxies arises at low redshifts, especially in the left two panels 
(z < 0.8). The AGNs at these redshifts show a displacement to 
warmer FIR colors compared to the general galaxy population, 
with a warmermedian value of the »Ji()o-i6o, and lack the long tail 
to cooler FIR colors that is seen among normal galaxies. In the 
two high redshift bins, the median values of AGNs and normal 
galaxies are more similar, though the AGNs still show a stronger 
wing on the warm side of the distribution. 

To quantitatively assess the level of AGN contamination in 
the FIR, we generate a suite of hybrid SEDs which combine 
AGN and SF templates picked to span the range of typical SED 
shapes found in real galaxies. We adopt three AGN SEDs from 
Mullaney et al. (2011) - a mean AGN SED, a typical high lu- 
minosity and a typical low luminos i ty SE D - as well as the 
mean QSO SED from Netzer et al.l (l2007b . suitable for AGNs 
with QSO-like luminosi ties. For the SF templ ates, we consider 
the five basis SEDs from'Mullane v et al.l(l201lh and add an M82- 
Uke and Arp220-like SED from flie SWIRE template library 



available from |http://sites.google.com/site/decompir/| 
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Fig. 3. Distributions of mioo-i60 for COSMOS PACS-detected AGNs (shaded histogram) compared to the general population of 
PACS detected galaxies in the COSMOS field. Each panel is a different redshift bin, indicated in the upper right corner of the panel. 
Also in the upper right is a measure of the KS probability Pks that the two distributions share a common parent distribution. The 
distributions are most different at z < 1 but become more similar at higher redshift. Note, systems with likely AGN contamination 
are not excluded from this plot. 



dPolletta et alj 120071) to incorporate strongly star-bursting and 
ULIRG-like systems. We create AGN-dominated hybrid SEDs 
by adding an AGN and a SF template scaled to give a fractional 
contribution of AGN Ught > 50% at a rest-frame wavelength of 
60 /vm. 

In each panel of Fig.|2]we plot a shaded polygon that encom- 
passes the range of observed-frame colors of these hybrid tem- 
plates after they have been appropriately redshifted to the upper 
and lower redshift bounds of each bin. The templates that are 
fully AGN dominated are typically very warm in both OT24-160 
and mio()-i6o, and many lie outside the plots. The cooler end of 
the polygons are populated by hybrids with ~ 50% AGN frac- 
tions. 

In all the panels, the polygons lie well away from the small 
black points which mark the location of the bulk of PACS de- 
tected galaxies. This is expected since most PACS detected 
galaxies are not likely to host an active nucleus, even those that 
may have obscured AGNs missed by X-ray surveys. However, it 
is also evident that only a few AGNs lie within these polygons 
at any redshift. By and large, most AGNs can safely be assumed 
to be dominated by cool dust emission in the FIR. 

The handful of AGNs that do lie in the AGN-dominated 
polygons are typically quite luminous. In addition, there is a 
weak but general trend for the more luminous AGNs to lie closer 
to the polygons in all the redshift bins. This implies that some 
small level of AGN contamination does enter into the FIR, even 
around and beyond 60 yum, and simply assuming that the rest- 
frame 60 ;um luminosity in AGNs is always attributable to star- 
formation is probably incorrect. It is clear, however, that the vast 
majority of PACS detected AGNs are dominated by SF-related 
emission in the FIR. 

Can a low level of AGN contamination in the FIR lead to 
the different wioo-ieo distributions between AGNs and inactive 
galaxies? Two features of Fig. |2] suggest otherwise. Firstly, note 
that the shaded area of strong AGN contamination is most offset 
from the the locus of the data points in the lowest redshift panel 
and closest to the locus in the highest redshift panel. This is sim- 
ply because the effects of AGN contamination on mioo-ieo are 
least pronounced at low redshifts, for reasons already discussed. 
Yet, it is in these low redshift bins where the most significant 
differences in color distributions are seen in Fig. |3] While the ex- 



clusion of AGN dominated objects does make the distributions 
more similar, substantial differences still remain, especially the 
lack of a cool wing to the AGN mioo-i60 distribution. Secondly, 
close examination of Fig.|2]reveals that objects that make up the 
spread to warmer FIR colors are not always luminous AGNs, 
as one would expect if contamination was the reason for the 
spread. In addition, the part of the color-color diagram occupied 
by 'warm' AGNs at low redshifts also contains some X-ray un- 
detected objects, implying that AGN activity is not a necessary 
prerequisite for having warm FIR colors. The lack of a close re- 
lationship between FIR color and AGN luminosity, coupled with 
an inverted trend with redshift, suggests that the offsets are not 
related to AGN activity, but are intrinsic to the population of 
galaxi es that host AGNs. 

In'Santini etal.1 (1201 2h . we show that X-ray AGNs at z < 1 
exhibit slightly enhanced mean SFRs compared to normal galax- 
ies of the same stellar mass, while the offset declines or disap- 
pears at higher redshifts. It is possible that the very processes 
that are responsible for the enhancement of mean SF activity in 
lower redshift AGNs may also affect their typical FIR colors. We 
discuss this further in §5. 

We turn now to AGNs that were not detected in both PACS 
bands. They account for the vast majority of X-ray selected 
AGNs (* 70% in the GOODS fields, more than 90% in the 
COSMOS field). PACS-undetected AGNs are less star-forming, 
and therefore, on average, more AGN dominated at a given AGN 
luminosity than the PACS-detected AGNs. We can study the FIR 
properties of PACS-undetected AGNs only through their mean 
stacked signal. In the lower panels of Fig. |2] we compare the 
mean colors of the PACS undetected AGNs to the PACS detected 
AGNs. The stacked points come from individual measurements 
in all three fields, while the black points here are only from the 
COSMOS dataset and serve only to place the stacked points in 
context. The stacks were restricted to objects that were detected 
in the 24 fim band, in order to allow an accurate estimate of the 
mean 24 fim flux. As before, we include the polygons that delin- 
eate regions in the diagrams occupied by AGN-dominated SEDs 
at 60 jum rest. 

By and large, the colors of the PACS stacked AGNs lie out- 
side the shaded polygons, implying that the stacked signal is 
not dominated at 60 yum by AGN light. However, as expected 
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they tend to lie closer to the polygons than the bulk of detected 
AGNs, clustering around the light edges where templates with 
AGN fractions of a few tens of percent would lie. Some or most 
of the stacks in any given redshift bin lie within the region of 
AGN dominance. In general, the scatter of stacked points lies 
closer to the polygons in the higher redshift bins, as expected if 
the shift to warmer colors is a result of AGN contamination. 

In addition, the stacks of the more luminous AGNs (red- 
der points) have consistently warmer colors. Most of the vari- 
ation with AGN luminosity happens through the m24-i6o color. 
Systematic changes in mioo-i60 are not clearly discernable and 
the typical FIR colors of the stacked AGNs are consistent with 
the typical colors of the detected AGNs. This highlights the fact 
that FIR colors are not greatly affected by increases in AGN 
luminosity, unlike the accretion sensitive MIR SED. We con- 
clude that, at lower redshifts (z < I ) AGN contamination in the 
FIR is minor in both PACS detected and stacked samples and 
very rarely gets above 50%. However, at higher redshifts out to 
z = 2.5, the AGN can substantially affect mioo-ieo and, therefore, 
our estimate of the rest-frame 60 jt/m luminosity. 

4. Results 

Using measurements of Leo, the mean rest-frame 60 yum lumi- 
nosity from all 3 fields combined (Appendix B) , we will now 
study the relationship between FIR emission and SMBH accre- 
tion to z = 2.5. We first look at mean trends that come from 
combining fluxes from detections and stacks in the PACS bands 
to get a single mean value for L^o in each bin in redshift and X- 
ray luminosity. In the following sections, we extrapolate Lx to 
a bolometiic AGN luminosity {Lacn) usi ng the following rela- 
tion, which is based on Eqn. 5 from Mai olino et al.l (l2007h and 
taking a ratio of 7 between th e luminosity at 5 100 A and Lagn 
(iNetzer & Trakhtenbrotl2"007l) . 



(logLx- 11.78) ,^„^^ 
log Lag/v = -^ijYi ^ ^-^^^ 



(1) 



where Lx (the 2-10 keV band X-ray luminosity) and Lacn are in 
erg s"'. 



4. 7 . Mean trends: Uo against Lacn 

In Fig.m L(,() is plotted against Lacn- The plotted Lacn is the 
median value for objects in all three fields in AGN luminosity 
bin, while error bars in Lacn show the interpercentile range in 
AGN luminosity containing 80% of the sample - typical errors 
in AGN luminosity are smaller than the error bars shown here. 
Different colors are used to represent measurements in the dif- 
ferent redshift intervals. All the mean measurements from the 
Herschel fields have been tabulated in Table [1] In the plot, we 
include as well data for the local sample of Swift BAT AGNs 
(black poi nts and line). I ncluded as well in the figure is the rela- 
tion from iNetzerl (|2009|) (dashed line) and colored solid curves 
which come from our fits in §4. 1 . 1 . T his plot may be compared 
directly to Fig. 6 of'Sh ao et al.l (l20I0l) and is analogous to Fig. 5 
inlMuflaney et al. (201^! 

We also include a shaded zone which shows the range in 
the plot occupied by pure AGNs (no star-formation) with the 
SEDs described in §3.1. These SEDs are expected to span the 
majority of true AGN SED shapes in nature. The tight coiTela- 
tions between Ly an d 12.3 fj.m IR luminosity among local AGNs 
dGandhi et alj|2009l) are used to link the bolometric AGN lumi- 
nosity (derived from Lx) and the FIR luminosity (through the 



AGN IR SEDs). While, at low AGN luminosities, torus emis- 
sion will not contribute measurably to the mean FIR luminosity, 
it plays an increasingly greater role at high AGN luminosities. 
Indeed, quite surprisingly, our Leo measurements for the most 
luminous AGNs at z > 1 are consistent with their mean SEDs 
having low or no star-formation. On the other hand, the average 
of AGN+host SEDs should be above the pure AGN SED if even 
some objects have star formation. This may indicate that our 
ado pted intrinsic AGN SEDs are too FIR-bright, as suggested 
bv .Mor&Netzeil (12012!) . 

At redshifts below z - 0.8, the variation of Leo with Lacn 
shows a consistent behavior At low AGN luminosities, Leo is 
weakly correlated or uncorrected with AGN luminosity, and 
approximately constant over two orders of magnitude in Lacn- 
This constant FIR luminosity, which we call L^y, is redshift de- 
pendent and rises by a factor of « 3 between z - and z = 0.65 
(the mean redshift of galaxies in the 0.5 < z < 0.8 bin). 

At higher AGN luminosities, a change is seen. Beyond a 
certain transition in Lacn, ^-60 begins to coiTelate strongly with 
Lacn- This is most evident in the local BAT sample and among 
AGNs in the 0.5 < z < 0.8 bins, which both span a sufficiently 
large range in Lacn to bring out this trend. The slope of the cor- 
relation is comparable in both redshift bins. The transition AGN 
luminosity, which we call L^^^, increases with redshift, by about 
an order of magnitude between z = and z = 0.8. 

While lower redshift AGNs show a particular relationship 
between Leo and Lacn, the situation appears to change at higher 
redshifts (z > 0.8). The correlation between Leo and Lacn is 
weak or absent. The overall trend is much flatter than for AGNs 
at lower redshifts, implying that beyond z ~ 1 the FIR emis- 
sion of the AGN hosts is quite independent of the accretion lu- 
minosity of the nucleus. L^^^ is not well characterised at these 
redshifts. 

Having said this, we note that these inferences come primar- 
ily from measurements of the highest luminosity AGNs in our 
sample, which are significantly lower than expected from an ex- 
trapolation of the low redshift correlation. This measurement in 
the 0.8 < z < 1.5 redshift bin is rather uncertain since it depends 
on a combination of only four objects and is subject to large sta- 
tistical uncertainties. The measurement in the 1.5 < z < 2.5 bin 
is more robust. We have individually examined the optical and 
IR images of these luminous AGNs, as well as the X-ray spec- 
tral fits, to ensure that they are not affected by incoiTect multi- 
wavelength associations or large eiTors. Therefore, we believe 
that the flattening of the L(,q-Lacn relationship among luminous 
AGNs at high redshifts is a real effect. 

The increase in the mean Leo of low luminosity AGNs con- 
tinues to z = 2.5. The change in L^g accelerates between 

z ~ 1 and z ~ 2, to reach LgQ~ 10"*^ erg s"' at 1.5 < z < 
2.5. This is somewhat at odds with the increase in the star- 
formation density of field galaxies with redshift, which typi- 
cally rises sharply to z '^ 1 and flattens out towards higher red- 
shift dHopkins & Beacoiiill2006h. However, AGN are typically 
hosted by massive galaxies (iBrusa et al.l 120091: IXue et al.]l2010l : 
iRosario et al.ll201 H) which have star- formation histories that are 
differe nt from the bulk of s tar-forming galaxies at these red- 
shifts dPerez-Gonzalez et al.i 2008). Star-formation 'downsizing' 
requires that such massive galaxies go through most of their star- 
formation around or earlier than z ~ 2 and contribute less to the 
integrated FIR luminosity density of the Universe at z < 1 than 
at higher redshifts. Therefore, we would expect AGN hosts to 
have a higher rate o f increase of L^g at 1 < z < 2 than at z < 1 
(ISantinietalJl2012h . 
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Fig. 4. Mean vLv(60/im) (Leo) vs. Lacn of X-ray selected AGNs in 5 different redshift bins from the local Universe to z = 2.5. The 
colored data points are combinations of mean measurements in 3 PEP fields: GOODS-N/S and COSMOS, while the black data 
points come from our analysis of the SWIFT BAT sample. The solid colored lines are functional fits to the m ean measurements, 
as described in §4.1.1. The dashed line is the correlation line shown by AGN-dominated systems in Netzej (1200 9). The shaded 
region corresponds to the approximate Icr range exhibited by empirical pure-AGN SEDs. At low redshifts, a strong change in the 
mean trend exists as a function of Lacn, which disappears at high redshifts. The mean Leo of low-luminosity X-ray AGNs increases 
monotonically with redshift, mirroring the increase in the mean SFR of massive galaxies across redshift. 



4.1 .1 . Fits to the Leo - Lacn trend 



We model the trend between log Leo and log Lacn in each of the 
five redshift bins as a combination of a flat line (a constant value 
L^q) and a straight line with a positive non-zero slope a. The 
constant value is determined by the mean SFR of low-luminosity 
AGN hosts and is unrelated to accretion activity. We capture the 
correlation between log Lacn and log Leo as a linear relationship, 
essentially a regression line. In total, three parameters define the 
fit: LgQ, a and the intercept of the regression line. The intercept 
also defines L^^^, which we take to be the value of Lacn where 
the regression line ordinate = Lgg, i.e, the AGN luminosity where 
SF activity that is correlated with BH growth starts to dominate 
the total SFR of the host galaxy. 

We fit the measured data using a non-linear least-squares fit- 
ting procedure (CURVEFIT in IDL). Results are tabulated in 
Table|2l In the two lower redshift bins where a good fit is possi- 
ble (local and 0.5 < z < 0.8) we arrive at consistent values for 
the slope of the regression line (a ^ 0.75 ± 0.23). Lg^ increases 
strongly with redshift, as expected. 



Among the two higher redshift bins, the transition is not 
well-defined. A best-fit slope of the regression line is lower 
(a ~ 0.5 + 0.23) and L^^^ is about 1.5 dex higher than at lower 
redshifts. However, the quality of the fit is rather poor (reduced 
X^ Z 3.5) in both these redshift bins and, in fact, a single straight 
line with a variable slope fits the data points just as well. 

A correlation between AGN luminosity (or mass accre- 
tion rate) and host galaxy SF rate has been reported by 
several studies, for local low-lum inosity Seyferts (Netzej 
l2009t iDiamond-Stanic & Rieke 201 1) and hig h redshift AGNs 
and O SOs dLutz et al.. .2008: .Hatziminaoglou et alj 120101: 
Bonfield eral.ll201lh . lNetzeil (120091) compared the SFRs of lo- 
cal SDSS Type II and LINER AGNs and find that they lie on 
the same relationship as luminous PG QSOs (Netz er et al. 20071) 
and the mm-bright QSOs at z ~ 2 dLutz et alJ2008h . The slope of 
this relationship is approximately 0.8, q uite con s istent with our 
own estimate of a. The relationship in iNetzeJ (l2009l) extends 
linearly to low values of AGN luminosity (Lacn< 10"*^ erg s"'), 
while we find a very de finite turnover t o a flat relationship below 
L^P^, most Ukely since iNetzen ( l2009l) specifically concentrated 
on AGN-dominated galaxies. 
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Table 1. Mean Rest-frame 60 /im luminosities in bins of Redshift and Hard-band X-ray Luminosity. 



Redshift bins 


Bins in De-absorbed 2-10 keV band X-ray luminosity Lx 






COSMOS 


41-42 


42-43 


43-44 


44-45 


45 


-45.5 


0.2-0.5 


43.17-44.04(18) 


43.72-43.93(100) 


44.07-44.33(34) 


— (2) 


- 


-(0) 


0.5-0.8 


-d) 


44.03-44.24(73) 


44.19-44.42(123) 


44.73-45.12(13) 


- 


-(0) 


0.8-1.5 


-(0) 


44.39 - 44.96 (28) 


44.59 - 44.68 (434) 


44.97-45.11 (134) 


44.57 - 


- 45.09 (4) 


1.5-2.5 


-(0) 


— (0) 


44.78-44.98(152) 


45.30-45.39(295) 


45.27 - 


45.62 (23) 


GOODS-S 


41-42 


42-43 


43-44 


44-45 


45 


-45.5 


0.2-0.5 


44.02-44.49(11) 


-d) 


— (1) 


-(0) 


- 


-(0) 


0.5-0.8 


43.90-44.14(30) 


44.04-44.45(19) 


43.94-44.35(11) 


-d) 


- 


-(0) 


0.8-1.5 


43.99-44.61 (18) 


44.40 - 44.93 (27) 


44.43-44.75(23) 


44.88-45.15(3) 


- 


-(0) 


1.5-2.5 


-(3) 


45.02 - 45.39 (20) 


44.66-45.05(32) 


44.87-45.13(5) 


- 


~(0) 


GOODS-N 


41-42 


42-43 


43-44 


44-45 


45 


-45.5 


0.2-0.5 


43.88-44.41 (12) 


42.95-43.81 (5) 


-(2) 


-(0) 


- 


-(0) 


0.5-0.8 


44.23-44.59(22) 


44.03-44.34(10) 


44.36 - 44.67 (5) 


-(0) 


— 


~(0) 


0.8-1.5 


44.53-44.77(33) 


44.69 - 44.98 (47) 


44.51-44.71 (39) 


44.28 - 44.75 (6) 


— 


-(0) 


1.5-2.5 


-(3) 


44.89-45.42(16) 


44.72-45.11 (44) 


44.08-45.27(5) 


- 


-(0) 



Notes. Mean Rest-frame 60 /jm luminosities are in units of log erg s 
the mean measurement in each redshift and X-ray luminosity bin. 



, Numbers in parentheses indicate the number of galaxies that contribute to 



Table 2. Functional fits to the Trends between Leo and Lacn in bins of Redshift 



Slope of Regression Line a left free 


Redshift bins 


L* ' 


a^ 


Lf 3 


x~ 


X (single fine fit) * 


Local 
0.5-0.8 
0.8-1.5 

1.5-2.5 


43.56 ±0.12 
44.18 ±0.08 
44.60 ± 0.09 
45.08 ± 0.08 


0.75 ±0.22 
0.84 ± 0.30 
0.58 ±0.18 
0.53 ±0.25 


44.23 ±0.13 
44.94 ± 0.33 
45.99 ± 0.64 
46.46 ± 0.82 


1.53 
0.41 
7.58 

7.57 


N/A 
N/A 
7.65 
6.97 




Slope 


of Regression 


Line a fixed at 0.78 




Redshift bins 




a 


1 '^ 

^AGN 


X' 





Local 


43.57 ± 0.08 


0.78 


44.27 ±0.11 


1.14 


0.2 - 0.5 


43.94 ± 0.09 


0.78 


44.82 ± 0.41 


6.47 


0.5-0.8 


44.17 ±0.07 


0.78 


44.89 ± 0.22 


2.42 


0.8-1.5 


44.62 ± 0.04 


0.78 


45.98 ±0.17 


5.61 


1.5-2.5 


45.12 ±0.06 


0.78 


46.51 ±0.21 


5.39 



Notes. See §4. 1. 1. *'' Mean Ljo of low-luminosity AGNs *^' Slope of Regression line valid for luminous AGNs ■'' Transition AGN luminosity: Ljo 
is significantly correlated with Lagn at z < 1 beyond this AGN luminosity '*' A single straight line fit between log Lagn and log Leo with a variable 
slope and normalisation. 



If we assume that the slope of the regression line is constant 
across redshift, as suggested by Netzer (2009), then the fits to 
our lower redshift data points gives the following relation: 



the relationship between SF and AGN activity between lumi- 
nous AGNs at low and high redshifts, which becomes important 
around z ~ 1 . 



log 



i^eo 



1044erg s- 



= log 



6.3 X 1044erg s' 



0.78 



/2) ''^■2.. Distributions of Uo against La 



We rerun our fits now fixing the regression line to Eqn. 2 and 
derive a second, more refined set of values for LL, which are 

' I I I I oO' 

also listed in Table |21 In Fig. 14] we include lines which show the 
expected trends in each of our five redshift intervals as solid col- 
ored lines. The data points for AGNs in the three lower redshift 
bins follow the trends quite well, since the trends were calibrated 
on these points. At high redshifts, the high Lagn measurements 
lie systematically off the expected trend, by as much as 3-4cr 
for AGNs at z ~ 1 and by around 2cr at z > 1.5. Leo for these 
highest luminosity AGNs appear to be comparable to, or per- 
haps slightly higher than, Lg^ at their corresponding redshifts. In 
other words, the mean FIR luminosity of the brightest AGNs at 
z > 1 are not very different from those of low-luminosity AGNs 
or normal massive galaxies at these redshifts. This may be con- 
trasted with the rather drastic change in the Leo with Lagn among 
AGNs at lower redshifts. This suggests that there is a change in 



Till now, we have concentrated on the mean FIR properties of 
the AGN population. Underlying these mean measurements is a 
distribution of FIR luminosities which can have considerable in- 
trinsic scatter We turn now to an examination of the Leo distribu- 
tion by considering separately properties of AGNs both detected 
and undetected in the PACS imaging. 

In Fig. |5] we plot vLy(60;um) against Lagn for PACS- 
detected AGNs on a diagram similar to that used for the mean 
measurements in Fig.|4] We also indicate the location in this di- 
agram of the stacked measurements for PACS undetected AGNs 
using colored numbers which list the number of objects that 
go into a given stack. Since the X-ray and PACS depths of the 
COSMOS field is considerably lower than the GOODS fields, 
we make separate panels in the Figure for COSMOS and the 
combined GOODS fields. 

If we concentrate first on PACS detected AGNs, a clear cor- 
relation can be seen between AGN luminosity and vLv(60/zm) in 
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Fig. 5. vLv(60/zm) vs. Lacw of X-ray selected AGNs in 4 redshift bins from z - 0.2 to z = 2.5. Galaxies in the GOODS and 
COSMOS fields have been separated into right and left panels respectively, as they differ greatly in both X-ray and PACS depths. 
Colored points are AGNs detected in both bands. Colored numbers show the location of stacked measurements and the number 
of objects that go into the stack. Dashed and solid lines, as well as shaded regions, are the same as in Fig |4] Note the apparent 
coiTelation in the colored points in the left panel, which is absent in the right panel. The correlation is primarily driven by selection 
effects inherent with wide and shallow surveys. 



the COSMOS dataset if we consider AGNs at all redshifts to- 
gether. While this correlation has a slope that matches that of the 
regression line that we derive above for high-luminosity AGNs, 
it is offset from the line by approximately an order of magni- 
tude in vLv(60;um). The coiTelation disappears, however, when 
one considers PACS detected AGNs in the GOODS fields, even 
though there is a large degree of overlap in Lagn between the 
two samples. The correlation among COSMOS AGN is not real 
- it is driven piimaiily by strong selection effects in the AGN 
sample with redshift. At high redshifts, the rare X-ray and FIR 
luminous AGNs which are not found in low redshift bins due to 
smaller survey volumes start to be represented significantly in 
the COSMOS sample. This effect can be seen even within a sin- 
gle large redshift bin (for e.g, the green points at 0.8 < z < 1.5) 
though the apparent coiTelation is weaker. The much smaller sur- 
vey volume of the combined GOODS fields ensures that lumi- 
nous AGNs are under-represented at all redshifts allowing for a 
more uniform AGN sample, which weakens the aforementioned 
selection effects and, by extension, the apparent correlation be- 
tween vLv(60yL(m) and L^gn- We suggest that relationships be- 
tween star-formation and nuclear activity in AGNs in the litera- 
ture that are based on individual detections from shallow, large- 
area surveys may be strongly aff ected by such selection effects 
(e.g.. lHatziminaoglou et alj|2010h . 

In general, the stacked measurements are up to an order 
of magnitude less luminous than the weakest PACS-detected 
sources (see also Fig. IB. lb . A substantial fraction of AGNs, 
especially at lower redshifts, are hosted by quenched galaxies 
(ISantini et al.ll2012 ). Since star-formation is very weak to negli- 
gible in such hosts, they are expected to lie well below the PACS 
detection limit at all redshifts and therefore contribute only to 
the stacked measurements. 



4.3. Mean trends: Nh against Lagn 

We turn now to the variation of L6() with X-ray obscuration, pa- 
rameterised by the column density of absorbing hydrogen Nh- 
The major merger model for AGN fueling and evolution predicts 
that most BH growth occurs when the AGN is enshrouded by a 
dense dust y obscuring envelop e at the center of a merger rem- 
nant (e.g.. ISanders et al]|1988h . If mergers are responsible for 
fueling most of the X-ray AGN population, as so me evolution- 
ary models propose (e.g., ISomerviUe et al.ll2008h . there should 
be definite trends between the obscuration of AGNs and their 
mean SFR. 

We explore this in Fig. |6] where we plot Leo against Nh in 
the four redshift bins that we used in earlier figures. To ensure 
a clean AGN sample, only sources with Lx> lO'*^ erg s"' were 
used to make this plot. At each redshift bin, we split the AGNs 
into six bins in obscuration: log Nh (cm^^) of < 20.5, 20.5-21, 
21-22, 22-23, 23-24, > 24. The plotted value of Nh for each 
data point is the median value for all AGNs in that bin and the 
error bars show the interpercentile range containing 80% of the 
sample. Measurements from all three fields have been combined 
in this Figure. 

In general, the trend between Lgo and Nh is weak. At 0.5 < 
z < 0.8, L(,Q appears to increase slowly with Nh by about 0.4 dex 
across 3 orders of magnitude in the obscuring column, though 
unobscured AGNs also show a higher mean Leo- The measure- 
ment for the intermediate bin in Nh is consistent with the L^q 
of the AGN population at this redshift, while bins with higher 
obscuration are slightly enhanced compared to this value. This 
trend appears to be driven more by selection effects in Lx-Nh 
space than by true correlations between L(,o and Nh (see Fig. [1] 
and associated discussion) and is not particularly significant . At 
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Fig. 6. L(,o vs. Nh in 4 different redshift bins in the redshift range 
0.2 < z < 2.5. The dashed Hnes show L^ in each redshift bin, 
taken from Fig. |4] All objects with Nh consistent with Galactic 
absorption contribute to the point at A^^ ^ 10^°'^ cm"^, irrespec- 
tive of the different Galactic absorption towards the three differ- 
ent fields. Leo is relatively flat across a large range in Nh with 
no clearly discernible trends, especially among well-measured 
points. As discussed in §4.3, covariances between Lx and Nh, 
driven by selection effects in the population of obscured AGNs, 
can affect trends in this figure, but the general results shown here 
should be relatively robust, since we combine measurements 
from fields of very different X-ray depths. 



0.8 < z < 1.5 and 1.5 < z < 2.5, the increase of Leo with Nh is 
weaker, consistent with a flat trend. 

At every redshift, a subset of AGNs have Nh < 10^°-^ cm"^, 
consistent with a very low or pure Galactic obscuration. In gen- 
eral, Leo for these objects is comparable to or higher than those 
of more obscured AGNs. 

We note here our AGN sample differs considerably from the 
luminous QSO samples for which a link between X-ray obscu- 
ration and SFR has been previo usly been suggested (Pa ge et alj 



ration and y K has been previo usly been suggested (Fa ge et al 
l200lLl2004trStevens et aliJOOSl) . The bulk of our AGNs are typi 
cally below the knee of the X-ray LF, while those studies concen- 
trated on special populations of X-ray absorbed optical broad- 
line QSOs. In addition, the typical AGN luminosity of our X-ray 
sample increases with redshift (Fig.[T}. At z < 1, most AGNs that 
contribute to Fig. |6]are of low and moderate luminosity, while 
at z > 1, they are increasingly luminous AGNs where merger- 
driven processes may be expected to be most relevant. 

Another simple, model-independent way to broadly separate 
obscured and unobscured AGNs is through their Hardness Ratio 
(HR), defined as HR ^ (H -S)/(H + S) where H are the photon 
counts in the hard band (2-10 keV) and S are the photon counts 
in the soft band (0.5-2 keV). If the typical SFRs of obscured 
AGNs were significantly enhanced, as suggested by models of 
AGN fueling by gas-rich major mergers, then we expect them 
to be more frequently detected by PACS (or other FIR surveys) 
compared to unobscured AGNs (e.g.. Page 2011). 

In Fig. |7l we plot HR against redshift for AGNs in the 
GOODS-S field. Despite the lower number statistics, we use 
this field for this plot because the PACS data is the deepest 
here and the detection fractions are not swayed by the statis- 
tics of the most FIR luminous systems, typically strong star- 
bursts. The PACS detection fraction s, estimated using Bayesian 
binomial statistics (ICameronll20Ilh . of AGNs above and below 



our fiducial boundary between unobscured and obscured AGNs 
(solid curve: Nh - 10^^ cm~^) are comparable in all four of 
our redshift bins. Similar results are found in the GOODS-N and 
COSMOS fields: any differences in the detection fractions be- 
tween obscured and unobscured AGNs are at a 1 cr level and not 
significant, given the typical uncertainties in X-ray counts. 
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Fig. 7. The observed hardness ratio (HR) of GOODS-S X-ray 
AGNs as a function of redshift. Large red points show PACS 
detected AGNs and small black points are PACS-undetected 
AGNs. Objects detected only in the hard(soft) bands have HR - 
1 and HR - -\ respectively. The solid black line shows the ex- 
pected HR for a power-law X-ray spectrum with a typical photon 
index of F = 1.9 and an obscuring Hydrogen column density of 
Nh - 10^^ cm"^. Dotted vertical lines mark the 4 redshift bins 
used in this study. The numbers in each bin give the Icr range in 
PACS detection fraction (in %) for obscured X-ray AGNs (i.e, 
those above the solid line - upper set of numbers) and unob- 
scured X-ray AGNs (i.e, those below the solid line - lower set of 
numbers). A quick examination will reveal that PACS detection 
fractions of obscured and unobscured AGNs are similar, with 
only insignificant differences in any given redshift bin. 



5. Discussion 

Using deep X-ray catalogs and high-quality FIR imaging from 
Herschel/PACS, we have estimated the mean FIR luminosities 
of AGNs across a wide span in redshift and SMBH accretion 
luminosity. We have studied the variation of the 60 /vm specific 
luminosity with Lagn and find that Leo is independent of Lacn 
at low accretion luminosities, but can show a strong correlation 
with Lagn at high accretion luminosities. The strength of the 
correlation depends on redshift - while strong at low redshifts 
(z < 0.8), the correlation appears to weaken or disappear at high 
redshifts. In this section, we discuss the implications of our re- 
sults on our knowledge of the process of AGN fueling, its evo- 
lution with redshift and its dependence on the nature of the host 
galaxy. 

The rest-frame FIR is widely regarded as the preferred tracer 
of star-formation in galaxies since it is particularly sensitive to 
cool dust emission which dominates the output of star-forming 
regions. However, as we show in §3, the tail of FIR emission 
from hot dust heated by an active nucleus can influence the FIR 
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Table 3. Mean Rest-frame 60 fim luminosities in bins of Redshift and X-ray Obscuring Column. 



Redshift bins 




Bins in Column Density of Obscuring Hydi 


rogen Nh (log cm -) 






COSMOS 


19-20.5 


20.5-21 


21-22 


22-23 


23-24 


24-25 


0.2-0.5 


43.78 - 44.07 (45) 


43.12-43.49(4) 


43.71-44.01(36) 


43.90-44.17(46) 


43.96 - 44.33 (4) 


-(0) 


0.5-0.8 


44.20-44.45(77) 


44.20 - 44.68 (5) 


44.06-44.31(38) 


44.02 - 44.37 (76) 


44.61-45.08(13) 


-(0) 


0.8-1.5 


44.74 - 44.89 (248) 


44.86-45.22(13) 


44.58 - 44.72 (88) 


44.61-44.76(212) 


44.71-45.00(37) 


~(2) 


1.5-2.5 


45.18-45.31 (289) 


-(2) 


45.02-45.33(35) 


45.09-45.31 (84) 


45.16-45.41 (58) 


45.13-45.49(3) 


GOODS-S 


19-20.5 


20.5-21 


21-22 


22-23 


23-24 


24-25 


0.2-0.5 


— (0) 


— (0) 


-(0) 


-(0) 


— (1) 


-(0) 


0.5-0.8 


43.68 - 44.69 (7) 


-(2) 


43.59-44.33(3) 


43.87-44.31(11) 


43.80 - 44.07 (7) 


~(1) 


0.8-1.5 


44.49 - 45.33 (8) 


-(2) 


43.96 - 44.47 (8) 


44.38-44.83(19) 


44.18-44.72(14) 


~(2) 


1.5-2.5 


44.82 - 45.40 (8) 


-(0) 


<44.59 (6) 


43.88-44.78(17) 


44.54 - 45.23 (23) 


<44.60 (3) 


GOODS-N 


19-20.5 


20.5-21 


21-22 


22-23 


23-24 


24-25 


0.2-0.5 


-(2) 


-(0) 


— (1) 


43.39-43.54(3) 


-(0) 


— (1) 


0.5-0.8 


44.06 - 44.29 (3) 


-d) 


— (2) 


44.22-44.51 (5) 


— (2) 


— (2) 


0.8-1.5 


44.74-45.12(19) 


44.13-44.44(3) 


44.04-45.29(11) 


44.40 - 44.64 (22) 


44.51-44.68(30) 


42.66 - 44.59 (7) 


1.5-2.5 


< 45.04 (11) 


-d) 


< 45.44 (4) 


44.73-45.44(14) 


44.82 - 45.23 (27) 


44.86-45.36(8) 



Notes. Mean Rest-frame 60 )im luminosities are in units of log erg s 
the mean measurement in each redshift and Nh bin. 



. Numbers in parentheses indicate the number of galaxies that contribute to 



emission in AGN host galaxies, and in some cases, may domi- 
nate the 60 fj.m luminosity of such systems. While the likelihood 
of AGN dominance in the FIR depends positively on the lumi- 
nosity of the AGN, this is by no means a simple one-to-one rela- 
tionship. Galaxies that are luminous in the FIR (PACS detected 
AGNs) are usually SF-dominated with a small fraction that are 
AGN-dominated. FIR weak galaxies (undetected in the R\CS 
imaging) are, on average, SF-dominated as well. 

5.1. Two regimes of AGN-galaxy co-evolution? 

Our results show that, among z < 1 galaxies with relatively lu- 
minous AGNs {Lacn<: lO'*'* erg s"'), the mean nuclear activity 
closely tracks the bulk of star-formation. At lower AGN lumi- 
nosities, on the other hand, star-formation is unrelated to AGN 
luminosity. This suggests two regimes in the relationship be- 
tween global star-formation in AGN hosts and nuclear activity. 

The first regime, applicable to low luminosity AGNs, re- 
quires a disconnect between instantaneous AGN accretion and 
host star-formation across the galaxy. The processes that drive 
and regulate most of the star formation in such host galaxies 
do not directly fuel the SMBH. The second regime, applicable 
to high luminosity AGNs, is one in which the SMBH accretion 
tracks the total SFR of the host galaxy. Here we expect direct 
coordination between global SF processes and the fueling of 
the AGN, over timescales that are short compared to the duty 
cycle of AGN activ ity (~ 10^"*^ yr. Martini & Weinberg 2001; 
iHopkins et all2005h . 

One of the best candidates for a process that directly links 
global SF and SMBH accretion are galaxy mergers, signatures 
of w hich are relatively common among local AGNs (iKoss et alJ 
l2010i) . Morphological studies of low -redshift star- forming QSOs 
also find freq uent m erger signatures ([Canalizo & Stocktonl200lt 
lUmitia et al]l2008l: IVeilleux et al.ll2009l) . OSO hosts frequentlv 
show bluer colors than ty pical old red elhptic als, in support of 
their recent merger origin (iJahnke et al.ll2004l) . However, a sub- 
stantial fraction of QSO hosts show characteristics of undis- 
turbed disks, inconsistent with a pure merger origin for al l QSOs 
(iBahcall et alJl997l:lDunloD et al.l2003l:lGuvon et al .'2006). The 
question of whether the most luminous AGNs are always asso- 
ciated with mergers is still somewhat open. Part of the reason 
for this uncertainty is that the structural identification of bona- 



fide mergers is not entirely straightforward (iLotz etal.l201 iL and 
references therein). Studies of moderate luminosity X-ray AGNs 
find no excess of merger signatures compared to inactive galax- 
ies (e.g. Grogin et al. 2005), and only weak trends wi th AGN lu- 
minosity (.Cisternas et al.,,201 U.Kocevski et alll2012h . However, 
these studies are usually Umited to AGN luminosities below 
^ACN' with very limited statistics at higher luminosities. In ad- 
dition, our results seem to suggest that the role of mergers as 
drivers of luminous AGNs may diminish at high redshifts, ex- 
actly where such studies have concentrated. 



5.2. The connection between SFR and Lacn in recent 
merger simulations 

It is instructive to test if the correlation observed among AGNs at 
z < 1, which we attribute to galaxy mergers, is indeed in agree- 
ment with models that trace the evolution of both star formation 
rate and accretion onto the black hole during a galaxy merger. 

We have approached this by re-visualizing results from a 
set of curre nt merger models. We use the models studied by 
IWuvts etaP (12010) which focus on progenitor galaxies with size 
and high gas content as expected for z~2 mergers, as well as 
a fully analogous set of models which enco mpass the larger 
sizes and lower gas content of local mergers (ICox e t al. 200q). 
In brief, all these models use the Gadget-2 smoothed particle hy- 
drodynamic code (Springel 2005) and prescriptions for star for- 
mation, supernova feedback, black hole growth and AGN feed- 
back to trace the evol ution of merging ga laxies. For more details, 
we refer the reader to lWuvts et al .1(1201 Ob and references therein. 
At z~2, 25 simulations cover a number of interaction geometries 
of compact high-z disks with gas fractions around 80%, while at 
z~0 30 simulations span different geometries for gas fractions 
5-20%. 

We caution that this suite of models samples only a small part 
of the possible parameter space of progenitor properties and in- 
teraction geometries, and that the treatment of star formation and 
accretion is approximate - the accretion rate follows the Bondi 
mechanism for the circumnclear conditions at the model's spa- 
tial resolution of order 100 pc. However, these models represent 
among the best attempts at this time to integrate SMBH accretion 
into galaxy mergers and can certainly allow a feasibility test of 
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Fig. 8. Average time spent by hydrodynamic merger models (see 
text) in different regions of the plane defined by AGN luminosity 
and star forming luminosity. Top: z~2 mergers of gas rich pro- 
genitors. Bottom: local mergers. The dashed line indicates the 
relation vLy(60f^m) k L°-^^ of Netzeret al. (2009). 



whether mergers produce an overall correlation between galaxy- 
wide SFR and Lagn- 

The two panels of Fig. [8] visualize the time that the mergers 
spend in different regions of the SFR/Lacn plane. These figures 
include tracks from all merger models together, separated only 
by their redshift regime. Any individual merger model generally 
takes a rather complicated path in the SFR/Lacn plane, suggest- 
ing that there is strictly no real tight synchronisation between 
the two processes. However, when taking the suite of models 
as a whole, it is clear that much tim e is indeed spent near the 
'correlation line' from iNetzen (l2009l) and §4.1, mostly when star 
formation and accretion decay after the merger Such behavior is 
seen among mergers appropriate for both low and high redshift 
progenitors. 

We refrain from making further detailed conclusions due to 
the limitations of this approach and since the scope of this dis- 
cussion is limited to a statement of consistency. While much fu- 
ture modelling work remains to be done to improve the tracing 
of mergers in the SFR/Lacn plane, results from these current 
models are consistent with mergers spending considerable time 
in the general area of the 'correlation line' of Fig. |4] 

5.3. Connections between global star-formation and AGN 
activity 

The existence of two regimes in the relationships between host 
SF and AGN activity suggests that there exists one or more 
key processes which modulate both star-formation and accretion 



onto the nuclear black hole. Since the two regimes are distin- 
guished by AGN luminosity, this implies that the main mode of 
AGN fueling changes between low and high luminosity AGNs. 
By this, we mean that the principal process (or family of pro- 
cesses) that drive the growth of AGNs is fundamentally different 
between the two regimes in AGN luminosity, and these two sorts 
of processes differ in their connection to global SF. 

Some insight into the processes that govern fueling at lower 
AGN luminosities comes from the well-studied population of lo- 
cal Seyfert galaxies, which account for most of the low and mod- 
erate luminosity AGNs in the local Universe (e.g., Simkin et al] 
[1980) . Seyfert host galaxies tend to be massive disks with sub- 
stantial s pheroids, though some fraction are also in more ellipti - 
cal hosts dKauffmann et al.ll2003Ll2007l:ISchawinski et alj|2010l) . 
While on-going major mergers are found among Seyferts, they 
account for only a minority of hosts. The statistics of galaxy 
morphology and color among local Seyferts strongly disfavors 
a major merger origin for the majority of local AGNs, since the 
big, stable disks in such galaxies could not have undergone a 
recent merger event. 

A plausible scenario, supported by such studies of the struc- 
ture and environment of local Seyfert galaxies is that SMBH 
fueling in low-luminosity AGNs is driven mostly by stochastic 
infall of gas from the circumnuclear region of the host galaxy 
^uld bears little connection to the extended SFR in the host 
dKauffmann & HeckmanI 120091) . Circumnuclear gas is replen- 
ished by several 'secular' processes that operate on the many 
galaxy-wide dynamical timescales. The large differences be- 
tween active periods of low-luminosity accretion and such long 
replenishment times disconnects the nature of gas in the galaxy 
at large and the more rapid. 

We find that the hosts of low luminosity AGN s are shghtly, 
but m easurably, different from inactive galaxies. ISantini et al.l 
(1201 2h find that they have a slightly enhanced mean SFR com- 
pared to normal galaxies of the same mass and in Fig. |3] we 
show that AGN hosts have warmer than typical FIR colors. This 
implies that, while the SFR in these hosts is not explicitly con- 
nected to AGN fueling, there is a preference for AGNs to lie 
in galaxies that are currently forming stars. The small shift to 
warmer colors among AGNs could be explained if AGNs are 
preferentially in star-forming galaxies that have more central or 
circumnuclear SF, where ambient gas densities and dust temper- 
atures may be a bit higher than in the outskirts of galaxies. It 
may be that the processes that trigger bursts of SF in the cen- 
ters of galaxies may also be related to gas inflow to the nucleus, 
though not in a synchronised manner This is supported by high- 
resolution si mulations of gas inflow to th e inner few parsecs in 
disk galaxies dHopkins & Ouataeral2010l) . 

Luminous AGNs {Lacn>Laqi^), on the other hand, are 

with 
20081 



fueled by rapid nuclear inflow of g as associate d 
major mergers of gas-rich galaxies (iMenci et al 



Volonteri et al."2003'; 'Hot?kins & BeacomI 120061: iHopkins etaP 
2008b; Somerville et al. 2008). Since they also inspire strong 
star-bursts, major mergers are a good candidate for a galaxy - 
wide process that directly correlates AGN accretion and host 
star- formation. In a sense, the merger-driven starburst-QSO con- 
nection may be thought of as a scaled up version of the SF- 
accretion relationship that may more commonly found among 
low-luminosity AGN, since in mergers the total SFR of the host 
galaxy is dominated by the nuclear starburst. 
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5.4. The role of major mergers at high redshift 

Turning now to z > 1, our results show that the correlation be- 
tween SF and AGN luminosity weakens or disappears at these 
redshifts. We proceed with an interpretation taking the points as 
they are measured, but caution that the interpretation is subject 
to the uncertainties associated with the measurement of Lf,o for 
the most luminous AGNs. 

Guided by our explanations in the above discussion of low 
redshift trends, we can develop a couple of explanations for the 
flattening of the high redshift relationship. It may be that the 
secular processes active in low-redshift Seyferts are capable of 
fueling most high luminosity AGNs as well. Several studies have 
established that gas-rich, turbulent disks are common among 
massive galaxies at high redshift, due to high inflow r ates o f 
cold gas (lErbetal.ll2006t iForster Schreiber et al.ll2006l l2009h . 
In such disks, a greater turbulent viscosity leads to a higher 
SMBH accretion rate (Boumaudet al. 201 lab), and possibly 
more frequent high Eddington-ratio bursts, without necessarily 
producing a direct connection between accretion rate and SFR. 
If greater turbulent accretion can fuel luminous AGNs, the role 
of mergers could be eclipsed at high redshift by fueling through 
such secular means. 

Having said this, current simulations still do not have the res- 
olution to link large stochastic bursts of accretion to enhanced 
secular inflow, so the idea developed above is not strongly mo- 
tivated by theory. If, instead, luminous AGNs at all redshifts are 
attributed to major mergers, then the lack of a correlation at high 
redshifts could occur if a regime of SF-AGN synchronisation, 
associated locally with the final coalescence of major gas-rich 
mergers, is lost among distant galaxies. This may happen if, for 
e.g, the interval of time between major mergers at high redshifts 
is short compared to the timescale of peak accretion, as cosmo- 
logical simulations seem to suggest. 

5.5. Implications for black hole-galaxy coevolution 

If the correlation between Leo and Lacn is indeed due to a dom- 
inant role of mergers among luminous AGNs, as supported by 
studies of local and low redshift populations, our results suggest 
a change in the relevance or nature of major galaxy mergers as 
drivers of AGN-galaxy co-evolution at high redshifts. We turn 
now to the implications of this in the context of galaxy evolu- 
tion. 

The most popular explanation for the form and tightness of 
the well-known local Mbh-o' relationship is that most black hole 
growth occurs in major mergers, where both the galaxy and the 
SMBH grow in lock-step. These high luminosity, high accretion- 
rate phases of SMBH growth are visible in their end stages as 
QSOs and evolutionary models of SMBH growth by mergers 
can explain simultaneously the form and evolution of empiri- 
cal SMBH mass functions, QSO luminosity functions, merger 
and red galax y density functions, clustering, and the X-ray back- 
ground (Hopkins et al.ll2008b a: Somerville et al. 2008) . 

Our results suggest that at low-z, most massive BH growth 
occurs in mergers or other possible processes that subject the 
host to violent relaxation. This maintains and tightens the mas- 
sive end of local SMBH scaling relations. 

At high-z, mergers may be less important. Many studies 
show that high-z galaxies are quite different from local ones. 
SF galaxies typically have large gas fractions an d high tur- 
bulent line-widths (Genzeletal. 2008; Forst er Sc hreiber et alJ 
120091: iDaddi et alJl2010l:lTacconi et al.ll2010l) . Such galaxies can 
have a high duty-cycle of unstable modes (clumps, bars, grand 



spirals) which can strongly torque gas. In addition, the larger 
magnitude of turbulent viscosity compared to local settled gas 
disks also plausibly leads to greater net 'secular' inflow to the 
nucleus. Simulations sug gest that sustained inflo w rates are en- 
hanced at these redshifts dBournaud et al.ll201 la!) . 

It is unclear if this larger net accretion rate to the nuclear 
regions translates into an increa se in the frequency of low- 
luminosity AGNs, as suggested by lBoumaud et al.l ( 1201 la ), or a 
boost across the board in the density of AGNs at all luminosities. 
The flattening of th e faint-end slope with redshift from X-ray lu- 
minosity functions ([Silverman et al.l l2008; Aird et al. 2010) sup- 
ports the latter explanation. If so, the flat L(,()-Lagn trend at 
z > 1 suggests that the role of mergers in fueling high luminos- 
ity AGNs may be replaced by global disk instabilities in gas-rich 
galaxies, which effectively removes the synchronisation between 
AGN accretion rate and SFR inherent to major mergers. Secular 
processes, which are important for low-luminosity AGNs at low 
redshifts, become dominant at all AGN luminosities at z > 1 . 

Regardless of the mechanism, our results suggest that the 
lock-step growth of stellar mass and galaxy properties becomes 
effectively decoupled from the growth of the black hole at higher 
redshifts. Most models of SMBH-galaxy co-evolution require 
that black holes grow in concert with their host galaxy (at least 
its spheroidal component). However, at the redshifts where both 
the stellar content and black holes in massive galaxies are put 
into place, the correlation between SMBH and galaxy growth 
seems to weaken or disappear This suggests that alternative 
ideas are needed to tie BH growth to galaxy growth, in order 
to yield the tight Mbh-o' relation seen today. Possible ideas are 
a strong role for feedback-regulated BH growth, coupled with a 
short lag between BH and bulge growth (Hopkins et al. 200S]) , 
or perhaps pure statistical e volution of both compon ents in a hi- 
erarchical merging context (IJahnke & Macci6ll201 ll) . 

6. Summary 

We have explored the variation of the mean far-IR luminos- 
ity of X-ray selected active galaxies with AGN luminosity and 
obscuration, using deep X-ray and Herschel/PACS datasets in 
three key extragalactic survey fields: GOODS-S, GOODS-N 
and COSMOS. Through the combination of both narrow -i-deep 
and wide+shallow surveys, redshift binning and the inclusion of 
stacked measurements for FIR-undetected subsamples, we avoid 
the effects of Malmquist bias which can produce incorrect corre- 
lations between SFR and AGN properties in FIR detected sam- 
ples taken only from wide-nshallow surveys. 

We gauge the importance of hot dust emission from the ac- 
tive nucleus at FIR wavelengths, using a suite of empirical AGN 
SEDs and SEDs of star-forming galaxies, chosen to cover most 
of the range of SED shapes expected for these two galaxy com- 
ponents. We compare the observed-frame MIR and FIR colors of 
AGNs and IR-bright inactive (X-ray undetected) galaxies in bins 
of redshift, finding that active galaxies generally span roughly 
the same color range as inactive galaxies at FIR wavelengths, 
though with a small tendency to host warmer FIR SEDs. We 
also demonstrate that a small minority of AGNs, typically the 
most luminous ones, are dominated by hot dust emission from 
the active nucleus. 

After appropriately combining measurements from all three 
fields, we find characteristic trends between AGN luminos- 
ity and FIR luminosity Leo. Low-luminosity AGNs display es- 
sentially no relationship between global star-formation rate, as 
traced by FIR luminosity, and the luminosity of the active nu- 
cleus. Instead, the Leo of such AGNs increases with redshift as 
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expected from studies of the SFR evolution of inactive massive 
galaxies. Black hole fueling in these systems is probably domi- 
nated by secular processes related to gas inflow and momentum 
transfer in galaxy disks and is unrelated to the overall growth of 
the host. 

On the other hand, luminous AGNs {Lacn> lO'*'*-^ erg s ') at 
low and moderate redshift (z < 1) show a strong correlation be- 
tween SFR and nuclear luminosity, implying a close relationship 
between the growth of the stellar component of the host galaxy 
and the growth of the SMBH. We interpret this as evidence for 
the importance of gas-rich major mergers in fueling high accre- 
tion phases in AGNs at these redshifts, consistent with the results 
of modern hydrodynamic simulations. 

At higher redshifts (z > 1), we find that the correlation 
flattens or disappears, implying a weaker connection between 
galaxy and black hole growth. This suggests that that the role 
of major mergers as drivers of black hole-galaxy co-evolution 
evolves with redshift. If so, this has important implications for 
co-evolutionary models and the evolution of the Mbh-c relation- 
ship. 

We find very little dependence between the mean FIR lumi- 
nosity and X-ray obscuring column Nh in X-ray selected AGNs, 
over five orders of magnitude in Nh- The population of obscured 
AGNs do not distinguish themselves in terms of their current 
SFR and various mechanisms are may be responsible for obscur- 
ing the nuclear X-ray emission of AGNs, with some certainly 
unrelated to global star-formation in the host galaxy. 
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Appendix A: FIR luminosities: detection, stacking 
and measurement 

As a direct tracer of the FIR emission we concentrate on the 
mean luminosity vLy{6QiJ.m), estimated at a rest-frame wave- 
length of 60 //m. This choice is a compromise between a wave- 
length long enough to avoid most of the AGN contamination (see 
discussion below) and at the same time short enough to be sam- 
pled by PACS 160 fim observations even at the highest redshifts 
considered in this work. 

To preclude any assumptions about the SED shape, we com- 
puted vLv(60/zm) though a log-linear interpolation of PACS 
fluxes, after convertin g them to luminosities following the ap- 
proach of IShao et all (gOlO). For the GOODS-S dataset, which 
includes 70 fiva observations as well, we interpolated between 
the two PACS bands bracketing rest-frame 60 /im. We checked 
that wider wavelength coverage in GOODS-S (the addition of 
70 i-im observations) does not introduce any bias in the analysis: 
when interpolating rest-frame 60 //m luminosity only from 100 
and 160 /^m PACS photometry, our results remain unchanged, 
and individual vLv(60;um) do not differ by more than 50%. 

We also explored fitting PACS fluxes with a typical IR tem- 
plate (e.g. from Chary & Elbaz 2001 or Dale & Helou 2002 li- 
braries) to derive vLy(60jjm). We found differences by a few tens 
per cent (only in 15% of the bins is the difference larger than 
50%), but the global picture is unaffected. Therefore, we adopt 
the simple interpolated vLy(60fim) estimate throughout, making 
no assumptions on the detailed SED shape. 

For sources fully detected (i.e. detected in both PACS bands 
used to interpolate 60 fim luminosity) we used the corresponding 



fluxes and computed individual vLJ,(60/im). For PACS sources 
that were completely undetected in both bands, we computed av- 
erage fluxes by stacking in a given Lx and redshift interval. We 
stacked at the X- ray positions on PACS residual maps using the 
iBethermin et al.l(l2010.) libraries. The use of residual maps, from 
which all detected objects were removed, avoids contamination 
by nearby brighter sources. PSF photometry was performed on 
the final stacked images. For each PACS band j, we then aver- 
aged these stacked fluxes with individual fluxes of partially de- 
tected objects (i.e. detected only in band j but not in the other 
one), weighting by the number of sources. These stacked and 
averaged fluxes in each band are used to get vL^^^'^'^ (6Qfj.m), in 
the same fashion as for the fully detected sources. 

The final 60 fim luminosity in each Lx and redshift interval is 
computed by averaging over the linear luminosities of detections 
and non-detections, weighted by the number of sources. Only 
bins with more than 3 sources are used in our analysis. 

Errors on the infrared luminosity were obtaine d by boot- 
strapping, in similar fa shion a s many curren t studies (IShao et al.l 
l20Tot iMullanev et all I20TI ISantini et all I20T2I) . A set of N 
sources, where N is equal to the number of sources per bin in 
Lx and redshift, is randomly chosen 100 times among detections 
and non-detections (allowing repetitions), and a vLy(6Qium) is 
computed per each iteration. The standard deviation of the ob- 
tained vLy{60fim) values gives the eiTor on the average 60 yum 
luminosity in each bin. The error bars thus account for both 
measurement errors and the error on the population distribution. 
They do not, however, account for cosmic variance. 

Appendix B: Combination across fields 

To reduce scatter and the effects of cosmic variance, we combine 
the estimates of the average vLy(6Qfim) in each bin in redshift 
and X-ray luminosity from all three fields into one mean number 
for each bin. 

The PACS maps in each of the fields differ in depth and area. 
Therefore, we begin by comparing the estimates of vLy{60iJ^m) in 
all three fields against redshift, as a simple test to ensure that the 
mean measurements in all three fields are in fact compatible. In 
the three upper panels of Fig. IB.li we plot the mean vLv(60/zm) 
from AGNs in each of the three fields separately. The points are 
colored by bins in X-ray luminosity. At a glance, one may gather 
that both the measurements and their redshift evolution in all 
three fields are very comparable. 

In the lower three panels of Fig. IB.li we plot separately 
the vLv(60//m) of PACS detected AGNs and those derived from 
stacks of the undetected AGNs. In addition, we have plotted the 
approximate rest-frame 60 fim luminosity limit as a function of 
redshift, derived from a log-linear interpolation of the 3cr PACS 
detection limits at 100 and 160 fim for each field. These lumi- 
nosity limits match reasonably well the lower envelope of the 
PACS detected AGNs in the Figure. The COSMOS detections 
are more luminous than those in the GOODS fields both at the 
luminosity limit, due to the significantly shallower PACS imag- 
ing, as well as at the upper end, because the larger area of the 
COSMOS field (5x larger than the GOODS fields combined) 
brings in more of the rare, very FIR luminous galaxies into the 
sample. The stacked points can be up to an order of magnitude 
fainter than the PACS luminosity limit. 

Due to the differences in the depth and noise properties of the 
PACS images in the three different fields, combining the data at 
the level of the maps poses difficulties. On the other hand, the 
measurements of FIR luminosity in each of the three fields show 
consistent values and behavior with redshift and Lx- Therefore 
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Fig.B.l. vLv(60yum) as a function of redshift for X-ray AGNs in the individual COSMOS, GOODS-S and GOODS-N fields. The 
points in all panels are colored by logarithmic bins in hard band X-ray luminosity (Lx), as shown. The upper panels show the mean 
Lgo estimated by combining fluxes from PACS detections and stacks. Note the consistency across all fields between the values and 
variation with redshift of the mean measurements. The lower panels show individual PACS detected AGNs and the stacked points 
for PACS undetected AGNs (in bins of redshift and Lx). The solid black lines show the approximate 3cr vLv(60/im) limit, derived 
from the flux limits of the 100 and 160 /im PEP photometry catalogs. 

we combined posf-/ac/o the estimates of the mean vLv(60;um) in 
each redshift and Lx bin. This was done by taking an average 
of the FIR luminosity in all three fields in each bin with a valid 
measurement, weighted by the inverse variance (1/eiT^) derived 
from the errors on the mean FIR luminosity. 
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